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Abstract 

During 2012-2013, a total of 4325 host-seeking adult ticks belonging to the genus Ixodes were collected from various 
localities of Hokkaido, the northernmost island of Japan. Tick lysates were subjected to real-time PCR assay to detect 
borrelial infection. The assay was designed for specific detection of the Relapsing fever spirochete Borrelia miyamotoi and 
for unspecific detection of Lyme disease-related spirochetes. Overall prevalence of 8. miyamotoi was 2% (71/3532) in Ixodes 
persulcatus, 4.3% (5/117) in Ixodes pavlovskyi and 0.1% (1/676) in Ixodes ovatus. The prevalence in /. persulcatus and /. 
pavlovskyi ticks were significantly higher than in /. ovatus. Co-infections with Lyme disease-related spirochetes were found in 
all of the tick species. During this investigation, we obtained 6 isolates of 6. miyamotoi from /. persulcatus and /. pavlovskyi by 
culture in BSK-M medium. Phylogenetic trees of 8. miyamotoi inferred from each of 3 housekeeping genes (glpQ, I6S rDNA, 
and flaB) demonstrated that the Hokkaido isolates were clustered with Russian 8. miyamotoi, but were distinguishable from 
North American and European 8. miyamotoi. A multilocus sequence analysis using 8 genes {clpA, dpX, nifS, pepX, pyrG, recG, 
rplB, and uvrA) suggested that all Japanese 8. miyamotoi isolates, including past isolates, were genetically clonal, although 
these were isolated from different tick and vertebrate sources. From these results, 8. m/yomofo/'-infected ticks are widely 
distributed throughout Hokkaido. Female /. persulcatus are responsible for most human tick-bites, thereby /. persulcatus is 
likely the most important vector of indigenous relapsing fever from tick bites in Hokkaido. 

Citation: Takano A, Toyomane K, Konnai S, Ohashi K, Nakao M, et al. (2014) Tick Surveillance for Relapsing Fever Spirochete Borrelia miyamotoi in Hokkaido, 
Japan. PLoS ONE 9(8): e104532. doi:10.1371/journal.pone.0104532 

Editor: Brian Stevenson, University of Kentucky College of Medicine, United States of America 
Received April 10, 2014; Accepted July 9, 2014; Published August 11, 2014 

Copyright: © 2014 Takano et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All accession number of sequences are 
deposited in DDBJ/EMBL/GenBank databases as the following accession numbers: AB900797-817 and AB904793. 

Funding: This work was supported in part by a grant for research on emerging and reemerging infectious diseases from the Japanese Ministry of Health, Labor, 
and Welfare (H24-Shinko-lppan-008) (http://www.mhlw.go.jp/seisakunitsuite/bunya/hokabunya/kenkyujigyou/hojokin-koubo-h25/gaiyo/19.html), and Program 
to Disseminate Tenure Tracking System, MEXT, Japan (http://www.jst.go.jp/tenure/index.html). The funders had no role in study design, data collection and 
analysis, decision to publish, or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 

* Email: a-takano@yamaguchi-u.ac.jp 

9 These authors contributed equally to this work. 



Introduction 

Borrelia miyamotoi, a member of the relapsing fever group (RF) 
borreliae, was first discovered from Ixodes persulcatus ticks and the 
rodent, Apodemus argenleus, in Hokkaido, the northernmost 
island of Japan [1]. Subsequently, B. miyamotoi has been found in 
Ixodes scapularis and Ixodes pacificus ticks in North America 
[2,3,4] and Ixodes ricinus in Europe [5,6]. Cases of human 
infections with B. miyamotoi were initially reported in Russia [7] . 
Following this report, several human cases have been confirmed as 
B. miyamotoi infections in North America [8,9] and Europe [10]. 
Recently, tick surveillance for B. miyamotoi was performed in 
Europe [11] and Russia [7]. The surveys showed that /. ricinus 
and /. persulcatus in the Eurasian continent consistently harbor B. 
miyamotoi with low prevalence. However, large-scale tick surveil- 



lance has not been conducted in Asian countries where /. 
persulcatus is distributed. 

In Japan, cases of Lyme disease (LD) are clustered in the 
northernmost island, Hokkaido, because its principal vector, /. 
persulcatus, is abundant throughout the island [12,13]. The tick is 
also a candidate vector for B. miyamotoi [1]. A previous study in 
Hokkaido detected B. miyamotoi from /. persulcatus ticks and from 
rodents [14]. However, widespread prevalence data of B. 
miyamotoi infections are still lacking for populations of ticks, 
including /. persulcatus and other genus Ixodes species. This basic 
information on prevalence of B. miyamotoi infections in ticks is 
urgently required for risk assessment of relapsing fever in 
Hokkaido. 

Ixodes pavlovskyi, a species closely related to /. persulcatus, is 
also distributed throughout Hokkaido [15,16,17]. Compared to 
the abundant /. persulcatus , I. pavlovskyi is a minor tick species. 
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Figure 1. Tick collection sites in this study. Ticks were collected in the areas shown by black dots. The gray shading shows district where ticks 

were collected during this investigation (see Table 1). 

doi:10.1371/journal.pone.0104532.g001 



Nevertheless, the close evolutionary relationship among /. 
persulcatus and /. pavlovskyi [16,17] raises concerns about their 
ability to maintain B. miyamotoi in nature. In this study, large-scale 
tick surveillance for B. miyamotoi was conducted in Hokkaido to 
estimate the infection rate of host-seeking adult Ixodes ticks. The 
tick-derived isolates of B. miyamotoi established in this surveillance 
were subjected to molecular analyses to characterize their genetic 
profile. The resultant field and laboratory data will serve as a 
baseline in understanding the epidemiology of B. miyamotoi in 
Japan. 

Materials and Methods 

Tick collection, DNA extraction and borrelial cultivation 

During the spring to summer seasons (April to July) of 2012 and 
2013, host-seeking adult ticks of/, persulcatus, I. pavlovskyi and /. 
ovatus were collected by flagging over vegetation from a total of 57 
forested areas in Hokkaido (Figure 1). In these areas, no specific 
permissions were required for collection of ticks, and this study did 
not involve endangered or protected species. The collection of/. 
ovatus was discontinued in 2013 because of its extremely low 
prevalence of B. miyamotoi infection. The number of ticks 
examined, and the prevalence of Borrelia spp. among them were 
calculated for each district. 

All of the ticks collected were subjected to a quantitative real- 
time PCR (qPCR) assay to detect spirochete DNA. The DNA for 



PGR templates was prepared from whole or half bodies of each 
tick by using sodium hydroxide (NaOH), as described previously 
[18]. Briefly, tick tissues including bacterial cells were lysed in 
50 uL of 25 mil NaOH for 10 minutes at 95°G. After adding 
4 (iL of Tris-HCL (1 M, pH7.5) for neutralization, the lysate was 
centrifuged at 4°C, and the resulting supernatant was used as 
template DNA for qPCR. A preliminary study suggested that the 
sensitivity of qPCR using the lysate was equal to that using DNA 
extracted by DNeasy Tissue Kit (Qiagen, CA, USA) (Data not 
shown). 

Parts of ticks collected were used to cultivate B. miyamotoi in 
modified Barbour-Stoenner-Kelly medium (BSK-M: using mini- 
mal essential medium alpha [BioWest, Germany] as a substitute 
for CMRL-1066) under microaerophilic conditions [19,20] (Table 
SI). The surface of ticks was sterilized by washing with 0.1% 
sodium hypochlorite, then with 80% ethyl alcohol. The washed 
tick was longitudinally bisected using a disposable knife (ELP 
No. 10, Akiyama Medical MFG. CO., LTD., Tokyo, Japan), and 
one half was inoculated into BSK-M medium. The remaining half 
was used to prepare the PCR template, as described above. Tick 
samples which were shown by qPCR to be positive for B. 
miyamotoi and negative for LD borreliae, were cultivated at 30°C 
for 4 weeks, and then the growth of spirochetes was examined by 
dark-field microscopy every two weeks. 
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Detection of borrelial DNAs from ticks 

As described previously [21], borrelial DNA in the tick lysates 
was detected by multiplex qPCR targeting the 16 S rRNA gene 
{16 S rDNA). Based on Barbour et al., a minor groove binder 
probe with an FAM-label (FAM probe) (Life Technologies 
Corporation, MD, USA) was designed for LD borreliae, while 
an minor groove binder probe with a VIC-label (VIC probe) was 
designed for RF borreliae, including B. miyamotoi. The qPCR was 
performed using Premix Ex Taq (Probe qPCR), (Takara Bio Inc., 
Shiga, Japan) according to the manufacturer's instructions. The 
qPCR was run on an ABI PRISM 7000 system, ABI StepOne 
system (Life Technologies Corporation), and LightCycler 480 
systemll (Roche Diagnostics, Basel, Switzerland). The PCR cycles 
were set at 40, and the reaction was performed in 12.5 uL (ABI 
PRISM 7000 system and StepOne system) or 10 U.L (LightCycler 
480 systemll) in single tubes or plate wells. The sensitivity and 
specificity of the qPCR was performed as described in the File S 1 . 

PCR, multilocus sequence analysis (MLSA), and 
phylogeny reconstruction 

Tick-derived isolates of B. miyamotoi were characterized by 
PCR-based DNA sequencing. After DNA extraction, using 
Wizard genomic DNA purification kit (Promega, WI, USA), 
Takara Ex Taq (Takara Bio Inc.) and KOD FX (TOYOBO Co., 
LTD., Osaka, Japan) DNA polymerases were used for all PCR 
amplifications, as recommended by the manufacturers. Target 
DNA fragments of the borrelial flagellin gene (flaB), glyceropho- 
sphoryl diester phosphodiesterase gene (glpQ), and 16 S rDNA 
were amplified as described previously [22] . MLSA was performed 
to further characterize the isolates of B. miyamotoi. For this 
purpose, the loci from 8 genes (clpA, clpX, nifS, pepX, pyrG, recG, 
rplB, and uvrA) were analyzed [23]. These primers were designed 
by genome assembling of 4 RF borreliae; Borrelia hermsii 
(CP000048), Borrelia duttonii (CP000976), Borrelia turicatae 
(CP000049), and Borrelia recurrentis (CP000993). The PCR 
thermal conditions to amplify the 8 genes were as follows: 95 °C for 



30 s, followed by 35 cycles of 95°C for 10 s, 50°C for 30 s, and 
72°C for 60 s. All PCR products were purified by High Pure PCR 
Product Purification Kit (Roche Diagnostics), then directly 
sequenced using an ABI Prism 3130 or 3130x1 Genetic Analyzer 
(Life Technologies Corporation). All of the PCR primers used in 
this study are listed in the Table S2. 

Raw sequence data were assembled by Sequencher 5. 1 (Gene 
Codes Corporation, MI, USA). Sequences of each gene were 
aligned by ClustalW (v. 1.6), and Neighbor-joining trees were 
generated by 1000 bootstrap repetitions under the Kimura 2- 
parameter by MEGA 5.2 [22,24]. All positions containing 
alignment gaps and missing nucleotides were eliminated only in 
pairwise sequence comparisons (the pairwise deletion option was 
used). 

Statistical tests for prevalence data 

Prevalence data of borrelial infections in tick groups from 
different categories (e.g. species, sex and locale) were analyzed by 
the Fisher's exact tests. A p-value less than 0.01 was considered to 
be statistically significant. 

Results 

Sensitivity and specificity of the qPCR 

The sensitivity and specificity of the qPCR were examined, as 
shown in the supporting information. The limit of the DNA 
detection consistendy observed was a minimum of 10 plasmid 
copies (data not shown). Regarding the specificity of qPCR, in the 
case of the FAM probe, 85% (11/13) of LD borreliae were 
detected, whereas Borrelia valaisiana Am501 T was barely 
detectable, Borrelia lusitaniae PotiB2 T was not detected, and 
none of the RF borreliae were detected. On the other hand, the 
VIC probe detected B. miyamotoi, B. hermsii, and Borrelia 
coriaceae Co53 T . However, for species with a 1 bp mismatch with 
the probe, such as B. duttonii Ly and Borrelia sp. AGRF [25], 
amplification was delayed, resulting in a reduced final concentra- 
tion that just barely reached our threshold. None of the reptile 




Figure 2. The prevalence of B. miyamotoi and LD borreliae in Ixodes ticks in Hokkaido. The asterisk shows significant difference: * P<0.01, 

by Fisher's exact test. 

doi:10.1371/journal.pone.0104532.g002 
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Figure 3. Phylogenetic analysis of RF borreliae based on 16 S rDNA of Borrelia spp. The phylogenetic tree of 16 S rDNA was constructed. 
The phylogenetic branches were supported in >70% by the bootstrap analysis. The bar indicates the percentage of sequence divergence. Sequences 
in this study were shown in bold type. If possible, clone or strain name, isolation source, and country were described in the case of B. miyamotoi. 
Numbers in parentheses indicate Accession Numbers in GenBank. 
doi:1 0.1 371 /journal.pone.01 04532.g003 



associated borreliae were detected by either probe. Therefore, the 
FAM probe widely detects B. burgdorferi sensu lato, which has 
been isolated from Ixodes ticks in Japan, including pathogenic 
borreliae of humans, such as B. garinii, Borrelia bavariensis and 
Borrelia afzelii, and non-pathogenic borreliae, Borrelia japonica, 
Borrelia tanuki and Borrelia turdi, whereas the VIC probe detects 
a narrow range of relapsing fever borreliae, of which the only one 
reported in Japan is B. miyamotoi. 

Prevalence of borreliae in Ixodes ticks 

A total of 4325 adult ticks were collected from various localities 
covering 11 out of the total 14 districts of Hokkaido (Figure 1). As 
shown in Table 1, ticks harboring B. miyamotoi were found from 8 
districts in northern, eastern and central Hokkaido. There was no 
significant statistical difference in the prevalence among the 8 
localities. Overall prevalence of B. miyamotoi in ticks of Hokkaido 



was 2% in /. persulcatus, 4.3% in /. pavlovskyi, and 0.1% in /. 
ovatus, respectively (Table 1). There was no statistical difference 
between /. persulcatus and /. pavlovskyi. The prevalences among 
both of these tick species were statistically higher than that of /. 
ovatus (Figure 2). The overall prevalence of LD borreliae was 
31.7% in /. persulcatus, 15.4% in /. pavlovskyi, and 42.6% in /. 
ovatus, respectively (Table 1). Co-infections of B. miyamotoi and 
LD borreliae in adult ticks were observed in 25 of 3532 /. 
persulcatus (0.7%), 1 of 117 /. pavlovskyi (0.9%), and 1 of 676 /. 
ovatus (0.1%). 

Phylogenetic analysis of isolated strains 

In this study, 6 isolates of B. miyamotoi derived from /. 
persulcatus and /. pavlovskyi were established in culture using 
BSK-M medium (Table 2). These isolates were utilized for 
molecular characterization, together with 5 isolates previously 
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Table 2. Borrelia miyamotoi strains used in this study. 



Strain name 


Isolation source 


Location 


Isolated year 


Reference 


HT31 T 


/. persulcatus, female 


Shiretoko, Hokkaido, Japan 


1990 


1 


HT24 


/. persulcatus, female 


Shiretoko, Hokkaido, Japan 


1990 


1 


FR64b 


Apodemus argenteus, whole blood 


Furano, Hokkaido, Japan 


1991 


1 


Hk004 


/. persulcatus, nymph 


Sibecha, Hokkaido, Japan 


1993 


1 


NB103-1 


/. persulcatus, nymph on Emberiza spodocephala 


Nemuro, Hokkaido, Japan 


1991 


1 


MYK1 


1. pavlovskyi, female 


Shibetsu, Hokkaido, Japan 


2012 


in this study 


MYK2 


/. persulcatus, female 


Nayoro, Hokkaido, Japan 


2012 


in this study 


MYK3 


/. persulcatus, female 


Obihiro, Hokkaido, Japan 


2012 


in this study 


MYK4 


/. persulcatus, female 


Simukapp, Hokkaido, Japan 


2012 


in this study 


MYK5 


/. persulcatus, female 


Takinoue, Hokkaido, Japan 


2012 


in this study 


Y13T1 


/. persulcatus, male 


Oumu, Hokkaido, Japan 


2013 


in this study 



doi:1 0.1 371 /journal.pone.01 04532.t002 



established in Hokkaido [1]. In the initial step of molecular 
analyses, DNA sequences of housekeeping genes (glpQ 16 S 
rDNA, and/Z<xB) were determined for all of the Hokkaido isolates, 
and there were no nucleotide substitutions in all 11 isolates. 
Therefore, the isolate HT31 (the type strain of B. miyamotoi) was 
selected as a representative for phylogenetic analysis. Neighbor- 
joining phylogenetic trees inferred from each of the housekeeping 
genes clearly showed that HT31 was clustered with B. miyamotoi 
from Russia, but was distinguishable from B. miyamotoi found in 
North American and B. miyamotoi detected from most /. ricinus 
ticks in Europe (Figures 3, SI, and S2). These three clusters of B. 
miyamotoi were temporally designated as "Siberian," "American" 
and "European" types, in this study (Figure 3). 

The molecular typing by MLSA using 8 genes (clpA, clpX, nifS, 
pepX, pyrG, recG, rplB, and uvrA) revealed that all of the 
Hokkaido isolates were identical to each other. Sequence similarity 
between the Hokkaido isolates (Siberian type) and American type 
of B. miyamotoi (LB-2001 [26], Acc. No. NC_022079) were 98% 
in concatenated sequences of the 8 genes (4776 bp), but the values 
ranged from 97.5 to 99.1% in individual genes. A phylogenetic 
tree constructed using concatenated sequences is shown in the 
supporting information (Figure S3). Nucleotide sequences of B. 
miyamotoi isolates from Hokkaido were deposited in DDBJ/ 
EMBL/GenBank databases as the following accession numbers: 
AB900797-807 and AB904793 for HT31 T , AB900808-16 for 
MYK1, and AB900817 for Y13T1. 

Discussion 

It is known that /. persulcatus principally serves as a tick vector 
for LD borreliae in Hokkaido. Hokkaido has an area of 
approximately 83,000 km 2 , and the total population of Hokkaido, 
according to the National Census 2010, is 5.5 million people 
(http://www.stat.go.jp/english/index.htm). However, 42% of the 
population inhabits the plains of Ishikari district (Figure 1). Similar 
to the taiga zone in Russia, /. persulcatus ticks are abundant 
mainly in forested areas of northern and eastern Hokkaido, which 
are thinly populated. Human tick bites and tick-borne LD, 
therefore, occur sporadically in Hokkaido [27]. In this study, we 
conducted large-scale tick surveillance for B. miyamotoi in 
Hokkaido, because human cases of B. miyamotoi infections have 
been confirmed in Russia [7]. Moreover, human B. miyamotoi 
infection, determined through retrospective analysis, has been 
reported in Hokkaido (http://www.nih.go.jp/niid/ja/relapsing- 



fever-m/relapsing-fever-iasrs/3877-pr4046.html). Our study dem- 
onstrated that B. miyamotoi is widely distributed throughout 
Hokkaido with low prevalence (0-3%) in /. persulcatus (Table 1, 
Figure 2). The prevalence is similar to those reported in Russian 
and European /. persulcatus [7,11]. Besides /. persulcatus, B. 
miyamotoi was found from 4.3% (within 0-6.1%) of/, pavlovskyi, 
and this is the first report of detection and isolation of B. miyamotoi 
from /. pavlovskyi. Moreover, we detected B. miyamotoi from an /. 
ovatus male. However, given that the prevalence in /. ovatus was 
significandy lower than that in /. persulcatus and /. pavlovskyi, the 
potential vectors of B. miyamotoi in Hokkaido are thought to be /. 
persulcatus and /. pavlovskyi. Female ticks of /. persulcatus are 
responsible for most human tick-bites in Hokkaido [27]. The 
abundance of B. miyamotoi-iaiecleA female ticks is, therefore, 
directly linked to human risk. Although the prevalence of B. 
miyamotoi is relatively low in rural communities of Hokkaido, B. 
miyamotoi is important for public health as an agent of fever and 
variable symptoms. Moreover, /. pavlovskyi can serve as an 
additional vector for B. miyamotoi. I. persulcatus and /. pavlovskyi 
are closely related morphologically, ecologically and phylogenet- 
ically [17,28]. The larvae and nymphs of both tick species have a 
wide range of feeding hosts, such as small mammals and birds. On 
the other hand, unlike /. persulcatus, adult /. pavlovskyi mainly 
infests birds. The potential for /. pavlovskyi to bite humans is still 
unclear, however, its public health significance is negligible given 
the far lower abundance. 

MLSA is an efficient tool to resolve phylogenetic relationships of 
LD borreliae on an intraspecies level [23]. In this study, we 
obtained for MLSA, 1 1 isolates, including 5 strains isolated in the 
early 1990's. These isolates were 100% identical in 8 loci, even 
from different tick and vertebrate sources. This result suggests that 
Japanese B. miyamotoi are clonal. The reason, however, remains 
unclear. Further analysis, such as the examination of MLSA using 
European B. miyamotoi or B. miyamotoi from the main island of 
Japan may contribute to the elucidation of its evolutionary history, 
e.g. clonal expansion. 

In previous studies, suspected human cases of co-infection with 
B. miyamotoi and LD borreliae have been reported from North 
America and Russia [6,7]. In this study, we also found co-infection 
of B. miyamotoi and LD borreliae in 3 tick species. Of 1118 /. 
persulcatus ticks which were infected with LD borreliae, 25 (2.2%) 
were co-infected with B. miyamotoi. It, therefore, seems likely that 
approximately 2% of LD patients might be co-infected with B. 
miyamotoi in Hokkaido. On the other hand, it is unclear when, 
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during feeding, B. miyamotoi is transmitted to vertebrates through 
saliva. Whereas LD borreliae are known to require several days 
(~48 hours) for transmission from ticks to vertebrates [29], soft 
tick-borne RF borreliae were efficiently transmitted through saliva 
within minutes [30]. In a previous study, we demonstrated that 
Borrelia sp. AGRF, which is phylogenetically related to B. 
miyamotoi, is maintained in the salivary gland of unfed hard-ticks, 
suggesting that the borrelial transmission to vertebrate hosts occurs 
rapidly when the ticks start to feed blood [25] . Thus, we speculate 
that B. miyamotoi may be transmitted to humans more rapidly 
than LD borreliae, thereby increasing the risk of infection. Further 
studies will be required to confirm this. 

In conclusion, the present surveillance showed that B. 
miyamotoi-iniected ticks are widely distributed throughout Hok- 
kaido, and that /. persulcatus is likely the most important vector of 
indigenous relapsing fever-like spirochetes transmitted by ticks in 
Hokkaido. Our subsequent genetic analyses also demonstrated 
that the Hokkaido isolates of B. miyamotoi are clonal. 

Supporting Information 

Figure SI Phylogenetic analysis of RF borreliae based 
on JlaB of Borrelia spp. The phylogenetic branches were 
supported in >70% by the bootstrap analysis. The bar indicates 
the percentage of sequence divergence. Sequences in this study are 
shown in bold type. The number in parentheses indicates 
Accession Number in GenBank. 
(TIF) 

Figure S2 Phylogenetic analysis of RF borreliae based 
on glpQ of Borrelia spp. The phylogenetic branches were 
supported in >70% by the bootstrap analysis. The bar indicates 
the percentage of sequence divergence. Sequences in this study are 
shown in bold type. The number in parentheses indicates 
Accession Number in GenBank. 
(TIF) 

Figure S3 Bayesian phylogenetic inference of concate- 
nated housekeeping gene sequences of RF borreliae. The 
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